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We present a theory which can explain the micromaser as well as its optical counterpart, the 
microlaser, for appropriate values of dissipative parameters. We show that, in both the the cases, 
the cavity radiation fields can have sub-Poissonian photon statistics. We further analyse if it is 
possible to attain a Fock state of the radiation field. The microlaser is precluded for such analysics 
due to the damping of its lasing levels making transitions at optical frequencies. Hence, we focus our 
attention on the micromaser and our exact simulation of the dynamics shows that it is not possible 
to generate a Fock state of the cavity radiation field. 



I. INTRODUCTION 

The subject of one-atom maser and its optical counter- 
part, the one-atom laser, has generated extensive interest 
after the recent experimental demonstrationa that they 
are capable of generating nonclassical states of the ra- 
diation fields [1-3]. In the one-atom maser experiment 
[1], two-level ^^Rb atoms in their upper Rydberg states 
of the transitions 63^p3/2 61^c?3/2(21.5065 GHz) and 
63^P3/2 61^^5/2(21.456 GHz) are pumped into a mi- 
crowave cavity at such a rate that, at most, one atom is 
present there at any time. So, there can be one of the 
two situations in the cavity: only one atom is present 
in the cavity or the cavity is empty of any atom. The 
cavity is tuned to one of the above two transitions in 
individual experiments and the sub-Poissonian nature of 
the cavity radiation field having variance less than that 
for a coherent state field has been inferred. It may be 
noted here that a well stabilized conventional laser can 
generate radiation fields that can at best be close to a 
coherent state field only. Sub-Poissonian radiation fields 
have photon distribution functions narrower than that 
for a coherent state field which has a Poissonian photon 
statistics. 

Whereas in the one-atom laser experiment [2], the 
pump rate of two-level 138Ba atoms in their upper states 
into a resonant optical cavity gives a stream of an aver- 
age number of atoms present there satisfying the con- 
dition (N) < 1.0. However, in order to have a sub- 
Poissonian radiation field, an uniform atom-field inter- 
action independent of cavity mode strucure is essential. 
The technique adopted by An el al [3] to improve their 
earlier setup [2] provides such a situation. Single-atom 
events in this cavity would be suitable for generating 
sub-Poissonian field in the optical regime since such ar- 
rangements in the microwave regime produced nonclas- 
sical fields [1] . 

In the following, we describe a unified formalism which 
shows that the one-atom maser (micromaser) as well as 
the one-atom laser (microlaser) cavity fields can have 
sub-Poissonian photon statistics. Our unified theory also 



explains the conventional laser dynamics. Further, the 
hope is that such sub-Possonian radiation fields would 
further shrink to bosonic Fock states. Our analysis below 
shows that this unified theory is not suitable for looking 
into these possibilities. Hence, we examine the possibil- 
ities of generating Fock states of the micromaser cavity 
field by an exact numerical simulation of the dynamics. 
In the case of microlaser, the damping of the atomic 
levels precludes such possibilities. 



II. THE FORMALISM 

We assume that atoms arrive individually at the cav- 
ity with an average interval i c = 1 /R where R is the flux 
rate of atoms. We have tc = t + tcav where r is the inter- 
action time, fixed for every atom, and tcav is the random 
time lapse between one atom leaving and successive atom 
entering the cavity, tc is the average of tc taken over a 
Poissonian distribution in time of incoming atoms. The 
cavity field evolves by this repititive dynamics from near 
vacuum as thermal photons in the optical cavity are al- 
most nonexitent. Thus, during r , we have to solve the 
equation of motion by taking into account the cavity dis- 
sipation as well as atomic damping, important at optical 
frequencies. Thus, we have 

p — —i[H, p] — K.[a)ap — 2apa) + pa) a) 

--f{S+S-p-2S-pS+ + pS+S-) (1) 

where H is the Jaynes-Cummings Hamiltonian [6] and k 
and 7 are the cavity-mode and atomic decay constants 
respectively, a is the photon annihilation operator and 
and S~ are the Pauli pseudo-spin operators for the 
two- level system. During tcav, the cavity field evolves 
under its own dynamics and is represented by 

p = — k(1 -I- nth){a^ap — 2apa) -f pa) a) 

—Kfithiao) P — 2a) pa + paa^) (2) 
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The method for obtaining a coarse-grained time deriva- 
tive for the photon number distribution P„ =< n\p\n > 
is given in detail in Ref. [4]. The steady-state photon 
statistics is then 

n 

Pn=Pol[ Vm (3) 

jn=l 

and Pq is obtained from the normahsation X]^o = ^■ 
The Vn is given by a continued fractions invoving the 
system parameters N = R/2k, the number of atoms 
passing through the cavity in a photon hfetime, n/g and 
7/5 where g is the atom-field coupling constant. Once P„ 
is obtained, we can get the laser intensity, proportional 
to (n), and its variance. This represents micromaser 
or microlaser photon statistics depending upon 
the relative values of the dissipative parameters 
and the coupling constatnt. We introduce another 
paprameter, the pump parameter D = ^/N gr useful for 
the description of laser characteristics. It is interesting 
to note that we can recover the conventional laser photon 
statistics from Eq. (2) by an appropriate choice of the 
above parameters. 



III. ONE-ATOM MASER DYNAMICS 

This involves the Rydberg levels of the active atoms 
having spontaneous lifetime tg = I/27 in the order of a 
fraction of a second. Whereas their flight time t through 
the cavity is of the order of microseconds. Hence we can 
safely put 7 = in the Eq. (1). We find from our for- 
malism that the normalized variance defined by 

V = V((n2) - (n)2)/(n) (4) 

is 0.5522 which is very close to the experimental find- 
ing V = 0.5477 in Ref. la near TV = 30 and values 
of the other parameters reported there. Thus we find 
that our formalism describes the micromaser dynamics 
quite accurately. It may be noted here that v = I for 
coherent state radiation fields. Thus, this novel device 
does generate radiation fields with sub-Poissoinian pho- 
ton statistics. 



IV. ONE-ATOM LASER DYNAMICS 

The lasing frequency here is in optical regime and, 
hence, the atomic lifetimes in these cases are short. 
Hence, we have to take 7 7^ for the microlaser dy- 
namics. For the typical parameters in the experiments 
in Ref. 3, we find that 7/5 = 0.1. Prom our photon 



statitics in Eq. (2), we obtain the laser intensity pro- 
portional to (n) as a function of the pump parameter 
D. The structures in (n) as D is varied for fixed N 
reflect the characteristics of the Jaynes-Cummings inter- 
action [6]. Soon after the threshold is attained at about 
D = 1.0, the photon number rises sharply. The reason is 
as follows. The field is almost in vacuum before the very 
first atom enters the cavity. Thus the atoms initially in 
their upper states contribute varying fractions of their 
energies to the cavity and at n = and D = it/2 the 
atoms get completely inverted. Thus (n) is peaked at 
about I? = 1.6 depending on k, 7 and N. For higher k 
and 7, the peak moves slightly towards higher D due to 
threshold being attained at higher D. We further find 
that < n >= around D =31.4, 62.8, 94.2,.... giving 
gr = IT, 27r, 37r,.... respectively. At such values of gir, the 
atom absorbs the photon it has emitted before leaving 
the cavity. 

The variance of the cavity field v increases sharply at 
about D = l.Q where (n) is also peaked. We find that 
near this value of D, the P„ is doubly peaked at n = 
and n N and this increases the variance in photon 
number. However, for slightly higher value of D, the 
cavity field is highly sub-Poissonian in nature. It also 
appears for further higher values of D. But, (n) is very 
small for such values of D due to increase in the inter- 
action time r which increases the infuence of atomic as 
well as cavity reservoirs on the atom-field interaction. It 
may be noted here that these characteristics are clearly 
different from that of conventional lasers. 



V. CAN THERE BE A FOCK STATE OF THE 
RADIATION FIELD? 

We sec that the one-atom dynamics is capable of nar- 
rowing the photon distribution function of the cavity ra- 
diation field compared to that for a coherent state field. 
This raises the question whether the distribution can be 
further reduced to a number state of the radiation field. 
The reason why such possibilities may arise is due to the 
fact that the reservoir influences are minimal in these 
dynamics. In fact, the analysis in Ref. [4] indicates this 
possibility in the micromaser dynamics when the condi- 
tion 

/(n) = -2nnth - 2N siv? {y/ngr) 

X exp[-7 - (2n - 1)k]t = (5) 

is satisfled. Clearly f{n) = only when nth = with 
gr = Tx 12. The later condition can be met by an adjust- 
ment of the atomic velocity with the help of the velocity 
selector in the micromaser setup [1]. But the other con- 
dition that the cavity temperature T = ^ K obviously 
cannot be satisfied. In the micromaser experiment [lb], 
the quality factore of the cavity is Q = 3.4 x 10^ and 
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its temperature is maintained at T = 0.3 K which gives 
thermal photons fith — 0.033. The Rydberg level in- 
volved in the atomic transitions have lifetimes of the 
order of a second which is much longer compared to its 
flight time of about r = 40 /xs through the cavity. These 
values of the parameters governing the dissipation dy- 
namics make their influences very small. In the case of 
one-atom laser, the lifetimes of atomic levels at optical 
frequencies arc short which makes the influence of atomic 
reservoir non- negligible. Hence, even though the thermal 
photons nth = in the optical cavity, the one-atom laser 
[2,3] would not be a right choice for this study. Thus, 
the micromaser turns out to be a suitable system to look 
into the possibilities of generating a number state of the 
cavity radiation field. 
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FIG. 1. Cavity photon distribution function at the exit of 
the 7000th atom. 



We then have to take into account the exact infuences 
of the reservoirs to acertain if a Fock state of the cavity 
field can actually be generated. Hence, instead of Eq. 
(I), we have to deal with 

p = —i[H, p] — k{\ + nth){a^ap — 2apa^ + pa^a) 

—Knth{ao)p — 2a)pa + paa)) (6) 

whenever a atom is present in the cavity. As mentioned 
earlier, a atom takes a time r to pass through the cavity. 
These atomic events are sepcrated by random durations, 
tcav^ during which the cavity evolves under its own dy- 
namics. Hence we set if = during tcav Processes like 
these stomic events seperated by random intervals are 
known as Poisson processes in literature encountered in 
various branches of physics, for example, radioactive ma- 
terials emitting alpha particles. A sequence of durations 
of such processes can be obtained from uniform deviates, 
also called random numbers, x generated using a com- 
puter such that < x < 1, and then by using the 
relation [7] 

tR = -p.\n{x) (7) 

where Ir = tcav + r. p, = 1/R where R is the flux rate 
of atoms. 

We have carried out numerical simulation of the dy- 
namics with the data taken from the experimental ar- 
rangements [1] in which g = 39 kHz and the t = 40 ps 
was one of the atom-field interaction times. In this case 
gr = 1.56 which produces Fock states of n photons where 
n satisfies sin gr^ n + 1 = in an ideal cavity {Q = oo). 
Since the experimantal arrangements are close to ideal 
situation, it was hoped that such Fock states could be 
attained experimentally. Indeed, such results have been 
reported in Ref. lb. However, our numerical simulations 
does not confirm these conclusions. Instead, it gives 
photon fields with very narrow distribution functions 
(sub-Poissoninan) centred about n. Figs. 1 and 2 dis- 
play distribution function P{n) narrowly centred about 
n = 14. 
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FIG. 2. P(n) vs n at the moment of the 9000th atom leav- 
ing the cavity. 



The reason for these results is simple. The cavity dis- 
sipation, although very small, effects the coherent atoim- 
field interaction and moreovc;r the randomness in tcav 
makes the photon distribution function fiuctuate all the 
time centred about n in addition to making it broader. 

In this experiment [1], the atoms coming out of the 
cavity are subjected to measuments from which state of 
the cavity field is inferred. The atoms enter the cavity in 
the upper \a) of the two states \a) and \b). The exiting 
atom is, in general, in a state 

IV) =a|a) + 6|6) (8) 

with Pa = |ap and pb = |6p arc the probabilities of the 
atom being in the states \a) and \b) respectively. Accord- 
ing to the Copenhagen interpretation of quantum me- 
chanics [8], this wave function collapses (or is projected) 
to either \a) or \b) the moment a measurement is made on 
it. Due to this inherent nature of quantum mechanics, a 
noise is associated with the measurement which is know 
as quantum projection noise [9]. Wc define the projection 
operator J = \a){a\. The variance in its measurement is 
given by 
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(AJ)2 = (j2)-(J)2=p^(l_pJ (9) 

We find that (A J)^ = only when Pa = 1 or 0. For the 
generation of a Fock state, it is necessary that the atom 
should leave the cavity unchanged in its upper state [4] . 
Hence, for such a situation we must have Po = 1 in which 
case (AJ)^ should be 0. We find from our numerical 
simulations that that pa = P{o.) is mostly about 0.8 
[Fig. 3] and, hence, (AJ)^ ^ always. This obviously 



indicates that the cavity field is in a linear superposition 
of Fock states giving a photon distribution function with 
the variance, defined in Eq. (3), i; > (For a Fock 
state V = 0). Indeed, we find that the v is about 0.5 in 
our calculations, presented in Fig. 4, indicating a sub- 
Poissonian nature of the cavity field. By itself, it carries 
a signature of quantum mechanics. 




We further notice in Fig. 4 that there are small fluc- 
tuations in V due to the fluctuations in P{n) [Figs 1 and 
2]. Also, V is nowhere near in Fig. 4. 
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FIG. 4. Fluctuations in v at the exit of successive atoms 
from the cavity. 



The steady-state photon statistics in Eq. (2) gives 
V = 0.597 which lies within the range of fluctuations in 
Fig. 4. This is due to the fact that the coarse-graining 



process involved in obtaining steady-state results ne- 
glects these small fluctuations (A detailed discussion can 
be found in Ref. [5]). Our exact numerical simulations, 
however, shows that they are crucial in the generation of 
a Fock state. We have carried out simulation until about 
10000 atoms passed through the cavity and carrying out 
the simulations any further would only be a repitition of 
the above fluctuations. 
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